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Abstract 

In this article we consider the frog model with drift on Z and inves¬ 
tigate the behaviour of the cloud of the frogs. In particular, we show 
that the speed of the minimum equals the speed of a single frog and 
prove some properties of the speed of the maximum. In addition, we 
show a limit theorem for the empirical distribution. 
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AMS 2000 subject classification: primary 60J10, 60K35; sec¬ 
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1 Description of the model and results 

We consider a system of interacting random walks, known as the frog model, 
on the one-dimensional lattice Z. In this model the particles performing the 
random walks are thought of as frogs. Initially, there is one active frog at 
the origin and one sleeping frog at every other site. The active frog performs 
a discrete time simple random walk jumping to the right with probability p 
and to the left with probability 1 — p where p > Whenever a sleeping 
frog is visited by an active frog, it is activated and starts a simple random 
walk itself with the same drift parameter p, independently of all other active 
frogs. It can wake up other sleeping frogs as well. 

Gantert and Schmidt prove in [6] that the frog model as described here is 
transient almost surely for every p > This means that the origin is visited 
only finitely many times by active frogs. Indeed, they consider a more general 
model with a random initial configuration of sleeping frogs and give criteria 
for recurrence and transience. In [T] and [2j Alves et al. prove a shape 
theorem in Z rf when the underlying random walk is symmetric. Transience 
and recurrence for frogs on Z* d is discussed in [16] . m and [5], for frogs on 
trees in [8] and [9], and in a more general setting in |T0] . Other aspects are 
discussed in in, m and [7]. An overview and a collection of problems up to 
the year 2003 can be found in [13] . 

In this article we investigate how the set of sites occupied by active frogs 
evolves over time. In particular, we consider the left and right front of this 
set and study their speed. In continuous time the right front is studied in 
[1] and [3] . 
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Let us introduce the model in a more formal way: Let (^)igz,fceN be an 
array of independent and identically distributed random variables such that 
P(X 1 ° = 1) = p = 1 - P(X° = -1). For i E Z, n E N 0 define S l n = Y2 =i X l 
The sequence (S'*,)neN 0 describes the trajectory of the frog initially at site i. 
It starts to follow that trajectory once it is woken up, which is a random 
time depending on the behaviour of the frogs activated earlier. We denote 
this activation time of frog i by A formal definition can be found in [Tj. 
The location Z l n of the frog initially at site i E Z at time n E No is given by 


ZL 


i for n < Ti , 

i + S l n _ T . for n > Tj. 


Let A n denote the set of active frogs at time n, i.e. A n = {i E Z: Tj < n}. 
Further, we dehne M n = ma Xj g ^ n Zf and m n = mirijg a u Z l n . Thus, M n 
describes the maximum and m n the minimum of the locations of the active 
frogs. We refer to M n and m n as the maximum and the minimum. One can 
show that there are non-zero constants n max and u m i n such that 

r M n 

Vmax = lim - a.S. 

n—¥ OO n 

r rn n 

^min — lim - a.S. 

n—>oo Tl 

The existence of u max is well known and is discussed in Section [2j and the 
existence of v m i n is part of Theorem 0 below. We call n max the speed of 
the maximum and w m i n the speed of the minimum. We study u max and u m i n 
as functions of the drift parameter p. First, we show that the speed of the 
minimum equals the speed of a single frog. 

Theorem 1.1. For p > \ the speed of the minimum exists and is given by 


Vm\n =2p-l. 


In the following two theorems we discuss some properties of the speed of the 
maximum. 


Theorem 1.2. The speed of the maximum is an increasing function in p. 
Theorem 1.3. For p < 1 it holds that u max < 1. 

In comparison to the last result note that for branching random walk on Z 
with binary branching the speed of the maximum equals 1 for every p>\, 
see e.g. m Theorem 18.3]. 

In addition to studying the behaviour of the minimum and the maximum 
we investigate the distribution of the active frogs. In the limit they are 
distributed uniformly inbetween the minimum and the maximum. To make 
this statement precise we rescale the positions of all active frogs at time n 
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roughly to the interval [0,1] and then consider the empirical distribution p n . 
which is defined for p < 1 by 


Vn{B) = 


| A 


E' 

i£.An 


(I'max-Vmin) 1 


F B ) 


for every Borel set B C [0,1]. Note that fi n is a random measure. 

Theorem 1.4. Almost surely, the empirical distribution p n converges weakly 
to the Lebesgue measure A on [0,1] as n —> oo. 


2 Proofs 


The existence of the speed of the maximum is proved using Liggett’s Subad¬ 
ditive Ergodic Theorem. Indeed, this theorem yields more information which 
we use throughout this article. We summarize it in the following Lemma. 


Lemma 2.1. For each p there exists a positive constant u max such that 


Furthermore, 


M n 

v ma , x = hm - a.s. 

n—>• oo Ti 


v = lim — = 

max ■ 

i —^oo l 


lim 

i—> oo 


E [Ti] 


= inf 

«SN 


E[Ti] 


a.s. 


(i) 


Proof. Let T t j denote the activation time of the frog at site j when initially 
there is one active frog at site i and one sleeping frog at every other site. An 
application of Liggett’s Subadditive Ergodic Theorem (see e.g. [12]) to the 
times (L;.j)ijez shows the existence of a positive constant u max such that ([!]) 
holds. For p = \ this is proved for a more general model by Alves et al. in |lj. 
In our setting their argument immediately applies to p > ^ as well. 

By a standard argument it now follows that lim n _ ) . 0O Mn exists almost surely: 
There exists a unique random sequence (fc n )neN with values in No such that 
Tk n <n < T kn+1 . Note that lirn^oo k n = oo. Hence, 


T n I )■ 

lim — 1 = hm —— 

n^oo n n^oo k n 


lim ] l 
n—>oo kn 


a.s. 


Obviously, k n — (n — T kn ) < M n < k n . This implies 


kn f kn \ -kin ^ k n 

n \ k n n J ~ n ~ n 


Taking limits yields the claim. 
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In order to prove Theorem [Ll] we compare the frogs initially on non-negative 
sites with independent random walks. The speed of the minimum of inde¬ 
pendent random walks can be computed explicitly which is done in the first 
of the following lemmas. Then it remains to deal with the frogs initially on 
negative sites. Luckily, they can be ignored due to the transience of the frog 
model proved in |6] by Gantert and Schmidt. 

We often need to talk about the frogs initially on negative sites. To keep the 
sentence structure simple we refer to them as the negative frogs. Analogously 
we speak of non-negative and positive frogs. Also for any k E Z the frog 
initially on site k is called frog k. 

Lemma 2.2. It holds that 

lim — min Si = 2p — 1 a.s. 

n-too n ie{-n,...,n} 

Proof. We only need to prove liminf n _ ) . 0O ^ min ie r n ^ n \ S l n > 2p — 1. For 
all £ > 0 we have 

P (— min S l n <2p — 1 — <2p—1 — e 

v ' 'i=—n 

/ s° 

< (2n + 1)P — < 2p- 1 -e 

\ n 

By Cramer’s Theorem the probability in the last term of this calculation 
decays exponentially fast in n, hence, it is summable. An application of the 
Borel-Cantelli Lemma and letting e —> 0 completes the proof. JIS. 


}) 


This result now enables us to prove a formula for the speed of the minimum 
of the non-negative frogs. 

Lemma 2.3. Let Af = {i > 0: Ti < n}. Then 


lim — min Z l n = 2p — 1 a.s. 

rwoo n i( zA+ 


Before proving Lemma 2.3 we make another short observation. Obviously 
Wiax is at least as big as the speed of a single frog, i.e. f max > 2p—l. In fact, 
this inequality is strict for all p > ^. For p = \ this fact is known from [I]. 

Lemma 2.4. For \ < p < 1 it holds that v max > 2p — 1. 

Proof. The key point in this proof is to notice that E[Xi] < E[T-f] holds, 
where Tf = inf{n € N: S® = 1} denotes the hitting time of the point 1 of a 


single simple random walk with drift 2p — 1. Hence, by Lemma 2.1 


= “ ^r 1 - E[ri1 < Era = 2F1 ■ 
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One can of course find better lower bounds for the speed of the maximum 
by estimating E[Tj] for i > 1, but this is not done in this article. 

Proof of Lemma \2.3\ It is enough to show liminfn^oc l n mim gj4 + Z l n > 2p—\ 
almost surely. 

In this proof we use a different but equivalent way of defining the movement 
of the frogs. For every j G Z define 



for n < Ti, 
for n > Tj. 


Note that (Zf) equals (Zf) in distribution. We now want to compare the tra¬ 
jectory (Z,\) ne pj 0 of each frog with the trajectory (S^) ne ^ 0 of the correspond¬ 
ing simple random walk. From time Tj onwards they move synchronously 
by the above definition. Therefore, we only need to compare their locations 
at time Tj. Note that Zf = i and define G = {i > 0: Sf, < i} to be the set 
of good frogs. Now i E Af n G implies S l n < Z l n for all n E N, i.e. all good 
frogs stay to the right of their corresponding random walk. Hence, 


min Z l n > min S l n 
i£Af 

> min S l n 


E ( S ‘n~ Z ‘n) 
ieG c nA+ 

E( s t-i). 

ieG c 


( 2 ) 


We claim that the set G c is finite almost surely. For p = 1 this is obviously 
true. For p < 1 it is enough to show that 


S l T . — i 

lim — f -— = 2p - 1 - n m ax a.s. (3) 

i—¥ oo 1 1 


since by Lemma 2.4 the last term is strictly negative and hence Sf. — i > 0 


can occur only for finitely many i > 0 almost surely. 

Note that (S , n ) n <x i is independent of the movement of the frogs up to time Ti. 
Therefore, Sf equals Sj,. in distribution. Using a standard large deviation 
estimate we get for every e > 0 


< 2p — 1 — = <2p-l-e S j < E[e _cTi ] < e~ ci 


where c = c(e,p) > 0 is a constant. By symmetry also > 2p — 1 + sj 

decays exponentially fast in i. An application of the Borel-Cantelli Lemma 
and letting e —» 0 thus shows 

St 

lim — 1 = 2p — 1 a.s. 

n—> oo Ti 
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Further, we know from Lemma 2.1 that lim^oo fr = v max almost surely. 
This proves equation (|3]) which implies that G c is finite almost surely. 
Therefore, the second term on the right side in inequality ([2]) is finite almost 
surely. Also note that it does not depend on n. Thus, 

1 1 

liminf — min Z l n > liminf — min S '* a.s. 

n^oo n i£ A+ n ^°° n i & A+ 


As A+ C { — n ,..., n} an application of Lemma |2.2| finishes the proof. 


Proof of Theorem 0 As shown in [5} Theorem 2.3] the frog model with 
drift, as considered here for p / is transient almost surely. This means 
that the origin is visited by only finitely many frogs almost surely. Therefore 
only finitely many negative frogs are ever activated. Hence, Theorem 1.1 


follows from Lemma 2.3 


Next we prove that the speed of the maximum is an increasing function in 
the drift parameter p. Though this statement might at first seem obvious, 
no direct coupling of the frog models for different drift parameters seems 
possible, since for smaller values of p more negative frogs will eventually 
be woken up, which might help in pushing the front forward. But we can 
ignore all these frogs without changing the speed of the maximum, similar 
to the proof of Theorem 0 This is shown in the next lemma. We therefore 
consider the frog model without negative frogs. It evolves in the same way as 
our usual frog model, but has another initial configuration. Here we assume 
that there is one sleeping frog at every positive integer, one active frog at 
0 and no frogs on negative sites. We denote the activation time of the i-th 
frog in the frog model without negative frogs by T~. 

Lemma 2.5. It holds that 


v 


-i 

max 


T+ 

lim —— a.s. 


n—>oo Tl 


Proof. We only need to prove limsup,^^ almost surely. First, 

we show that the speed of the maximum of all negative frogs in the usual 
frog model equals 2p— 1 almost surely, i.e. setting A~ = {i < 0: Tj < n} we 
prove that 

lim — max Z'l = 2p — 1 a.s. (4) 

n^oo n i( zA~ 

For p > ^ only finitely many negative frogs will ever be activated almost 
surely as remarked in the proof of Theorem 0 and proved in |6]. In this 
case equation Q is thus obvious. 

If p = then by symmetry the claim follows from Lemma 
Let E be the set of all positive frogs which are activated by negative frogs, 
meaning that at the time of their activation at least one negative frog is 


2.3 
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present. Since u max > 2p — 1 as proved in Lemma 2.4 and by equation Q 
the set E is finite almost surely. 

Hence, T = sup ig £ Ti ) is an almost surely finite random variable. For 

all i £ E we thus have T) + < Ti + T. Actually, this inequality is true for all 
j G No, which immediately implies the claim of the lemma. 

The inequality can e.g. be proven inductively. For i = 0 the inequality is 
obviously true as Tq = To = 0. Now assume that i G N and T+ < Tj +T 
holds for all 0 < j < i — 1. If i G E, there is nothing to show. Otherwise, let 
0 < k < i — lbe the (random) frog that activates the frog i in the normal 
version of the model. Then we have 


T+ < T+ + (Ti - T k ) <Ti + T. 


Note here that in both models the frogs follow the same paths, they might 
just be activated at different times. Ji§. 


Proof of Theorem 1.2. Using a standard coupling of the random variables 
(xi) ieZ,fceN we can achieve that X (p) is monotone decreasing in p. As 
»ma x(p) = hm n _ > . 0O ffr almost surely by Lemma 2.5 we conclude that v max (p) 


is increasing m p. 


In order to bound the speed of the maximum from above we prove an upper 
bound for the number of frogs in the maximum. We do this for a slightly 
modified frog model: Each time the maximum moves to the left we put a 
sleeping frog at the site that has just been left by the maximum. Hence, 
in this new model there is one sleeping frog at every site to the right of 
the maximum at any time. Further notice that, except at time 0, there are 
always at least two frogs in the maximum. We use the same notation as in 
the usual frog model, but add an index “mod” when referring to the modified 
model. Further, let a n denote the number of frogs in the maximum in the 
modified frog model. 

Lemma 2.6. For \ < p < 1 and all n € N it holds that 


E[a n ] < 


(2 ~ p)p 
(\-p)(2p-iy 


Proof. We prove bounds not only for the number of frogs in the maximum, 
but for every other site as well. Therefore, let a n (k ) be the number of frogs 
at location M“ lod — 2k for k,n E No- We prove by induction over n that for 
all o 

s (i-r)(2r-V ' (5) 

For n = 0 and n = 1 one easily checks that the claim is true. Assume that 
the claim holds for some integer tjGN. 
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First we show inequality (|5]) for k = 0. Distinguishing whether all a n particles 
in the maximum at time n move to the left or not in the next step one 
calculates 


E[a n+ i] = E [(1 - p) an (a„ + pa n ( 1))] 


+ E 


(i - (i ~pY 


pa n 


+ 1 


.1 - (1 -p) a 

= E [(1 - p) an (a n + pa n ( 1) - l) + pa n + l]. 


Note here that the expectation of a binomial random variable with para¬ 
meters p > 0 and k G N conditioned on being at least 1 is given by ■ 

Using a n > 2 yields 

E[a n+ i] < (1 - p) 2 E[a n + pa n (l) - l] + pE[a n ] + 1. (6) 

Inserting the induction hypothesis ([5]) in ([6]) the claim follows after a straight¬ 
forward calculation. 

For k = 1 an analogous calculation yields 


E[a n+ i(l)] = E[(l -p) an (pa n (2) + (1 - p)a n { 1))] 

pa n 


+ E 


(i-(i -pY 


dn. 


l - (l -p) a 
= E[(l -p) an \pa n { 2) - (2 p- l)a n (l) - a n )] 
+ E[(1 - p)a n + pa n ( 1)]. 


+ pa n 


(7) 


For k > 2 one gets 


E[a n+1 (k)} = E[(l — p) an (pa n (k + 1) + (1 - p)a n (k))] 

+ E[(1 - (1 - p) an )(pa n (k) + (1 — p)a n (k - 1))] 

= E[(l — p)“ n (pa n (k + 1) - (2 p - 1 )a n (k) - (1 - p)a n (k - 1))] 
+ E[(l — p)a n (k — 1) + pa n (k)\. (8) 

Thus, for k > 1 equations Q and ([8| imply 

E[a n +i(A:)] < p(l - p) 2 E[a n {k + 1)] + pE[a n (k)\ + (1 - p)E[a n {k - 1)]. (9) 

As before, inserting the induction hypothesis (|5]) into inequality Q com¬ 
pletes the proof. iHH 

Proof of Theorem Consider the event that in the modihed frog model at 
time n all a n frogs sitting in the maximum move to the left. Using Jensen’s 
inequality and Lemma |2.6[ we conclude that the probability of this event is 
bounded from below by 

r 1 nr 1 ( 2 -p)p 

E[(l — p) an ] > (1 - P ) E[an 1 > (1 - p) (l-r)(2p-i) . 
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Therefore, for all n G N 


E[I££? - T” od ] > 1 + 2E[(1 - p ) an ] 

(2-p)p 

> 1 + 2(1 - p) (i-p)Pp-i). 


Clearly, in the modified model, frogs are activated no later than in the normal 
version of the frog model. Thus, 


n 


E [T n ] > E[T“ od ] = \ E[T fc mod - Tf™?] > (1 + 2(1 - p) d-ripp-D ) n . 


fc=i 


(2-p)p 


By Lemma 2.1 we conclude 

E[T„ 


v mL = inf 

raeN 


n 


(2—p)p 

> 1 + 2(1 — p) ( 1 _ p)( 2 p _1 ) > 1. 


It remains to prove Theorem 1.4 The idea of the proof is quite simple: From 
the point of view of the minimum the front moves with a positive speed, but 
all the frogs only fluctuate around their locations with y/n, so basically they 
stay where they are. 

First, we show that for large enough times n all active frogs do not deviate 
much from their expected locations. To formalize this statement we define 
the set G n = {i E A n : | Z l n — E[Z^]| < n 3 / 4 }. 


Lemma 2.7. Almost surely, G n = A n for all n large enough. 


Proof. As A n C {—rz,..., n} we have 

P {An / G n ) = p( U (I Z l n~ nz l n ]\ > U 3/4 }) 

i(zA n 

n 

< E H\ z% n-n z i}\ >n 3/4 ) 

i=—n 
n n 

= E E P d Z n-E[4]|>n 3/4 |^ = A:)-P(T i = A:). (10) 

i=—n k =0 


Further, for every i G Z and 0 < k < n it holds that 

P(|Z; - E[Z;]| > n, 3 / 4 1Tj = k) = P(|^_ fc - EK_ fc ]| 


< 2 exp 

< 2 exp 


n 3 / 2 . 
4(n — k)J 
n 1 / 2 \ 

T"/' 


> n 3 / 4 ) 
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In the first inequality in the above estimate we use Hoffding’s inequality. 
Thus, (10) implies 


n l/ 2 nn nl/2 


G n ) < 2exp(^-— j ^2 = k) < 2(2n + l)exp 

i=—nk=0 

which is summable. An application of the Borel-Cantelli Lemma completes 
the proof. 

For e > 0 and x G [0,1] define 

J{i G Z: - (u max - e)n < * < ((2x - l)u max - e)n\ for p = L 
L n (x,£) = > 


{* G Z: 0 < i < (xw max - e)n} 


for p > 1 


and 


j {/ G Z: ((2x - l)v max + e)n<i< (w max - e)n} for p = 
R n {x,e) = < , 

I {* G Z : (rv max + e)n<i< (u max - e)n) for p > 5 . 


Lemma 2.8. For n Zarpe enough, i G L n (x,e ) D G n implies 


Z~ n V m j, | H 

(%as Wiin)^ 

whereas i G R n (x,e ) D G n implies 


< x, 


( 11 ) 


•^n Wiin^ . 

> X. 


('^max Wiin)^ 


( 12 ) 


Proof. For p = 1 note that by symmetry o m i n = —u max . Thus, (11) holds if 


and only if Z l n < (2x —l)u max ?r. Assume i G L n (x,e)nG n . A straightforward 
calculation shows 


Z l n < E [Z l n ] + n 3 / 4 = i + n 3/4 < (2® - 1) 




for n big enough. Analogously, one shows (12) in this case. 

For p > \ the proof works essentially in the same way as for p = \, but the 
estimation of E [Z\\ is less trivial. We have E [Z l n ] = i + (n — E[Tj])u m i n . For 
i G L n (x, e) D G n we thus get 


Zf <E[Z*]+n 3 / 4 = u min n + 




E [Ti] 


^min^max 


+ n 


3/4 


Lemma 


2.1 


yields that > iuf,; e N = f m ax- Hence, for n big enough 


Z n A V m i n Tl T (Wiax Wiin) T Tl 


3/4 


Wiin^ T ®(^max Wiin)?L 


10 















as claimed in (11). On the other hand, i E R n {x, e)nG n analogously implies 


P ^miiW T 


E[Ti] 


min 17 max 


— n 


3/4 


Since lim^oo = v^x and i tends to infinity whenever n does by the 
definition of R n (x,e), we know that + 5e for n big enough and 

a small constant 5. Therefore, 

Z n P U m i n Tl H (fmax — ^min £<5^min^max) Tl ^ . 

%ax 

Using i > (xu max + e)n and choosing 5 small enough finishes the proof. 


Proof of Theorem We need to show that lim n _ 5 . 00 /r n ([0, x]) = A([0,x]) 
for every x E [0,1] almost surely. 

Take a realisation of the frog model such that A n = G n holds for suffi¬ 
ciently large n, that lim, woo ^ = u max and lim^oo ^ = u min , and finally 
that A n D Z _ is finite. This happens with probability 1 as we have seen in 
Lemma 2.7 Lemma |2.1[ Theorem |1.1| and previous discussions about the 


transience of the frog model. Now fix x E [0,1] and e > 0 small. Lemma 2.8 
yields that, for n large enough, 


Pn([0 ,x]) > 


IA 


1 ir> n T (n. - n l L «( X ’ e )| 

I G n n T n \X, £) | — 


\A r 


n 


(13) 


For the last equation we used that L n (x,e ) C A n for sufficiently large n as 
lirrpi^oc = U m ax- The definition of L n (x,e) implies 


\L n (x,e)\ > 


j 2(xv max — e)n for p=\, 
\(xv max -s)n for p>\. 


-nr r = l v mlx for p = 5 , respectively lim^oo rfU = v 


-i 

max 


Further, linin-^oo ppr — j^max 1U1 f ~ 2 ’ lc °F c '- l > iVC1 .y iiiu inco ppy 
for p > \. Thus, the limit inferior of the last term in (13) as n —> oo is 
bounded from below by x — Since e > 0 was chosen arbitrarily we 

conclude 

liminf fi n ([0, xl) > x. 

71—> OO 

On the other hand, Lemma |2 . 8 | shows that, for n large enough, 


1 


Hn([0,x]) < y—-\A n \ (G n D R n (x,s)) \ = 1 - 


n \R n (x,e)\ 


\A n 


n 


(14) 


since A n = G n and R n (x,e) C A n for n big enough. By the definition of 
R n {x, e) we have 


\Rn(x,e)\ > 


2((1 — x)v max — e)n for p = \, 
((1 - x)v max - 2e)n for p > \. 
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Analogous to the above estimation this yields that the limit superior of the 
right hand side of (14) is bounded from above by x + As before we 

get, since e > 0 is arbitrary, 


limsup / u n ([0, x]) < x, 

n—i-oo 

which finishes the proof. 


3 Open Problem 

Simulations suggest that the speed of the maximum is a concave function in 
the drift parameter p. 

1 

0.9 

0.8 

M n 
n 

0.7 

0.6 

0.5 

0.5 0.6 0.7 0.8 0.9 1 

P 

Figure 1: Simulation of for n = 100000 

A heuristic argument might be the following. We expect that the number of 
frogs in the maximum converges to a stationary distribution a p for p < 1. 
Therefore, the speed of the maximum should equal 

£M*0(1-2(1 -?)*)• 

fee N 

If a p was independent of p. this would be a concave function. However, we 
believe that the dependence on p does not destroy the concavity. 

Conjecture 3.1. The speed of the maximum is a concave function in p. 
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